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We calculate, using a generalized parton distribution based formalism, the nonlinearity of the 
Rosenbluth plots that arise from hard two-photon exchange corrections that are not included in the 
classic calculations of the radiative corrections to electron-proton elastic scattering. 



I. INTRODUCTION 

Previous investigations HH2H3HH El El have shown 
that a possible explanation for the discrepancy in the 
measurement of the ratio of electric to magnetic proton 
form factor, measured through the Rosenbluth separa- 
tion technique [7. 8 9 \ and polarization transfer experi- 
ments 11T01ITT1IT21 , is due to a two-photon exchange pro- 
cess. Such process gives significant correction to the re- 
sult of Rosenbluth technique, while it gives a relatively 
small correction to the polarization transfer method. 

The usual strategy in analyzing experimental e-p elas- 
tic scattering data is to remove the classic radiative cor- 
rections 1 13 14i [T5l[T6] , most often using the Mo and Tsai 
expressions [15], and present the data in a form to be di- 
rectly compared to the lowest order expressions. For the 
differential cross section, one has the Rosenbluth expres- 
sion 

d<7 B orn 

dVL lab ~ e(l + t) dCl lah ' K ' 

where the "no structure" cross section c^g is well- 
known, t = Q 2 /(4M 2 ), and e = (1 + 2(1 + 
t) tan 2 9/2) , with M being the proton mass and 8 the 
electron lab scattering angle. The leading order reduced 
cross section is given in terms of the magnetic and elec- 
tric form factors, 

C7|°™ = G 2 1 (Q 2 ) + ^G 2 (Q 2 ). (2) 

The reduced cross section at fixed Q 2 is linear in e, and 
allows a separation of the two form factors. 

Additional, hard two-photon, corrections can change 
the normalization and e slope of the reduced cross sec- 
tion, and also introduce terms that are not linear in £. 
At high Q 2 , the G\ contribution to the reduced cross 
section is small, so that e dependent corrections to the 
big terms have a significant impact on the extraction 
of Ge/Gm- Calculations of these effects underlie the 
recent works flj 13 IH 13 |6| that at least partly rec- 
oncile the Rosenbluth and polarization transfer mea- 
surements of Ge/Gm- The same calculations also give 
non-linearity in e to the Rosenbluth plot. Experimen- 
tally 1 17, 20 1, the upper bounds on the non-linearity are 
tightening, and it has been noted that the bounds are 
now in the vicinity of what is expected from theoretical 
calculations. One can, as has been done [17], estimate 
the nonlinearity induced by the two-photon corrections 
from published plots. However, we wish here to pre- 



cisely quantify the non-linearity directly from the calcu- 
lation of the hard two-photon corrections, in particular, 
quantify it in two-photon calculations where the inter- 
mediate hadronic state is treated at a quark level J3J 2], 
described using generalized parton distributions. 

Recently, Chen et ah (21] have also considered non- 
linearity in the Rosenbluth plot. Chen et al. parameter- 
ize the two-photon contributions to the reduced cross 
section using functional forms that satisfy requirements 
based on crossing symmetry [22]. The values of their pa- 
rameters they obtain by fitting directly to data, and they 
show that the two-photon corrections can significantly 
change the e-slope of the reduced cross section while in- 
troducing only a modest nonlinearity. We here are doing 
a more specific calculation using the GPD based model, 
and get explicit predictions as to what the nonlinearity 
should be, in the context of this model. We can also ana- 
lyze the sensitivity of the predictions to different model 
GPDs and to different analytic representations of Gf or 
Gm- 

In the next section, after a few remarks about how 
two-photon exchange was treated farther in the past, we 
will record how the formulas for the reduced cross sec- 
tion are modified by the hard two-photon process, and 
then continue with a description of how to specify quan- 
titatively the non-linearity. In the third and last section 
we will show calculated results for the nonlinearity pa- 
rameter, and make some closing remarks. 




FIG. 1: Two-photon exchange contribution 



II. RADIATIVE CORRECTIONS AND FORMULAS 

Most radiative corrections to electron-proton elastic 
scattering were calculated long ago iT3llT4l|T5llT6I . No- 
tably not calculated in the far past were the full two- 
photon exchange contributions, Fig. [T] Since energy 
transferred by one photon can be returned by the other, 
it is possible for both photons to be quite energetic and 



hence probe proton structure quite deeply. The pro- 
ton structure knowledge needed to complete the two- 
photon calculation was not available then but is avail- 
able now, at least to the extent of having sufficient 
knowledge of proton form factors and / or generalized 
parton distributions that one can and should attempt a 
full calculation. At least two groups have explicitly done 
so EH5J. 

The two-photon corrections were not completely ne- 
glected in the past. If one or the other photon is 
soft, there will be infrared (IR) divergences, and these 
IR divergences must be calculated and included be- 
cause they cancel other IR divergences coming from 
bremsstrahlung corrections. Calculations of the two- 
photon corrections in the classic papers were hence done 
with approximations to the numerators of the expres- 
sions that are exact when one photon or the other is soft, 
but in general amounted to neglecting in the numerator 
terms of order of the loop momentum-squared. Thus, 
although there was an explicitly expressed hope that the 
neglected contributions would be small [14J, there is the 
real possibility that they could be large. 

We turn to a summary of the modern two-photon re- 
sults following [4J . Let cr be the measured cross section 
after the classic Mo-Tsai corrections have been applied. 
With the additional two-photon corrections, 
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+ ^G 2 E (Q 2 )^(l + mx) + <T RMrd . (3) 



The hard corrections are 
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where A, B, C, and Cr,/^ depend upon both Q 2 and e, 
and the characteristic integrals are 
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Each of the three GPDs that enter, W, W, and W, are 
evaluated at zero skewedness, for example, 



W = Hl{x,Z,Q 2 ) = W{x,0,Q 2 ) . 



(6) 



We picture the two photons interacting with a single 
quark, as in Fig. [2] The Mandelstam variables for the 
overall process are s, u, and Q 2 , and the evaluation is 
facilitated by using a frame wherein (external momenta 
are labeled in Fig. El 

\ 

V = n (P +P) = (P + >P~>P±) 




FIG. 2: Partonic scattering process of direct and crossed box 
diagram 



3+,(M 2 + Q 2 /4)/p+,0 ± 

(k' + k) = (rjp+,Q z /(4 VV +),Q x ^ 
-p = k-k' = (0,0, if±), (7) 



and one calculates 
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The momentum fraction x in the integrals for A, B, and 
C is for the active quark in the proton, and the Mandel- 
stam variabes for the electron-quark subprocess are s, u, 
and Q 2 with, 
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The elementary electron-quark amplitudes are obtained 
from electron-muon amplitudes 1 18 1, and their real parts 
are [4 ] 
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Given a GPD, one can now calculate the two-photon 
corrections. However, the GPDs are not at this time 
definitively known. We evaluate with two model GPDs, 
named the "gaussian" and the "modified Regge" Ifl9l , 
whose parameters are quoted in Ref . [4] . Both are con- 
strained to give the quark distribution functions in the 
appropriate limit, and fit as well as they can the form 
factors when appropriately integrated. The modified 
Regge GPD gives a better fit to the form factors, so 
would be the current best choice if one chose just one 
GPD. 

The curvature of the reduced cross section is seen in 
Fig. [5J which shows 
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FIG. 3: Fractional correction S = (a R - (7-| ™) /a R torn for two 
model GPDs and Q 2 = 6 GeV 2 . 



for Q 2 = 6 GeV 2 . Since a partonic calculation is valid at 
high momentum transfer, we require —u> M 2 , which 
leads to the low-£ cutoff on our curves. 

To quantify the nonlinearity, we, following [17J, fit 
the reduced cross section to a quadratic polynomial in 
e. With the Q 2 dependence tacit, 



o-R(e) 



Pi 



+ P2 



(12) 



where the p; are constant in £ (though functions of Q 2 ). 
If the curve is not quadratic, the p, will depend on the 
interval fitted. We will, at any Q 2 , fit to the region £j ow 
to 1, where £i ow is fixed by the boundary — u = M 2 , or 



M 1 



M 2 + Q 2 /2 



(13) 



Since a theoretical calculation yields a continuous 
curve, we can fit using a three-point gaussian method, 
the same fit method that underlies two-point gaussian 
integration. By way of reminder, if ones has a function 
f(x) defined for the interval — 1 < x < 1, one can ex- 
pand the function up to the second Legendre polyno- 
mial, 



f(x) = c Po(x) + ciPi(x) + c 2 P 2 (x) • 



(14) 



Noting that Pi is zero for x = X\ = and that P 2 is zero 
for x — x±2 = ±1/ V3/ one can solve for the q in terms 
of f(x) evaluated at three points, 
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(15) 



As a side note, the gaussian integral is the integral over 
the Legendre polynomial fit, 



dx f(x) = lag . 



(16) 



For a fit interval a < £ < b, one maps the above pro- 
cedure linearly using 



b + a b — a 
£ = — 1 — x or x 
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2 2 
The expansion maps to 

cr R (e) = a P (x) + a x Px{x) + a 2 P 2 (x) , (18) 
where x is given by Eq. ( [17) . The coefficients become 
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where 
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Simple algebra relates the coefficients p, for the expan- 
sion in (e — 1/2) of Eq. 1 12 1 to the q just given, 
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c 2 . 



We use = £i ow and b = 1, and have a fit procedure that 
is quite easy to code. 



III. CURVATURE PARAMETER RESULTS AND 
CONCLUSIONS 



The curvature coefficients P2 defined in Eq. (12) that 
follow from the two GPDs mentioned previously are 
shown in Fig. H] 

The Figure shows several form factor combinations. 
In addition to the Kelly form factors 1231 , we there are 
also plotted two other curves for each GPD. The second 
uses the magnetic form factor from Kelly and the electric 
form factor from the experimenters's linear fit 1 11 ], 



g e = (1-0.13 (q 2 - 0.04) 



(22) 



for Q 2 in GeV 2 . One sees that this change in Gg has 
little effect. Using instead the Brash et al. empirical fit 
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Nonlinearlity parameter, calculated with I u I > M 2 
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FIG. 4: Non-linearity parameter for c R <x 1 + pi(e— 1/2) + 
p 2 {£ - 1/2) 2 , fitted to calculated a R , with -u > M 2 . 



GPD dependence is quite noticeable. 

The survey of experimental Rosenbluth data by 
Tvaskis et al. |T7| obtains (p 2 ) = 0.019 ± 0.027 and a 95% 
confidence level upper limit |p 2 | max — 0.064. These are 
averages over p 2 obtained from many values of Q 2 , and 
nearly all of the individual p 2 (Q 2 ) have much larger un- 
certainty limits. The most notable individual points are 
from the Q 2 = 2.64 and 3.2 GeV 2 Jefferson Lab measure- 
ments which by themselves yield (p 2 ) = 0.013 ± 0.033, 
albeit these were measurements where the proton rather 
than the electron was observed 1 25 1 . 

Our nonlinearity coefficients calculated using the 
modified Regge GPDs are negative, and about 1 .4 stan- 
dard deviations away from the latter points. One should 
note again that the GPDs are not yet definitively known, 
and that the gaussian GPD, yields a significantly larger 
nonlinearity. Still, it seems clear that the experiments are 
at the point of providing a serious test of the two-photon 
exchange calculations, and in this regard we note that 
there is a dedicated Rosenbluth nonlinearity test run- 
ning in Hall C at Jefferson Lab. We look forward to the 
analyzed results. 



to G^vf [24 1 gives results barely distinguishable from the 
previous two curves; we think the same will be true of 
any fitted Gm that represents the data well. For some- 
thing more extreme, the third curve for each GPD uses 
the dipole fit for Gm supplemented by the above linear 
fit for Ge/Gm- Thus the Gm dependence of the Rosen- 
bluth curvature is slight for an up-to-date Gm fit, but the 
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